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The low-temperature physics of quantum many-body systems is largely governed by the structure 
of their ground states. Minimizing the energy of local interactions, ground states often reflect strong 
properties of locality such as the area law for entanglement entropy and the exponential decay of 
correlations between spatially separated observables. Here, we present a novel characterization of 
quantum states, which we call ‘local reversibility’. It characterizes the type of operations that are 
needed to reverse the action of a general disturbance on the state. We prove that unique ground 
states of gapped local Hamiltonian are locally reversible. This way, we identify new universal features 
of many-body ground states, which cannot be derived from the aforementioned properties. We use 
local reversibility to distinguish between states enjoying microscopic and macroscopic quantum 
phenomena. To demonstrate the potential of our approach, we prove specific properties of ground 
states, which are relevant both to critical and non-critical theories. 


I. INTRODUCTION 

Gapped ground states define quantum phases of matter at zero temperature. Even though they occupy a 
tiny fraction of the possible many-body Hilbert space, these states manifest a rich and diverse structure. Standard 
examples are states with local order-parameter such as paramagnetic and ferromagnetic ground states, the superfluid 
and insulator ones in bosonic and fermionic many-body systems, etc. Other instances, such as quantum Hall and 
quantum spin liquids, can arise because of more subtle orders that can be established in the system. A central 
goal of condensed matter theory is to understand their structure and how it relates to the physics of different 
phases [1, 2]. A natural approach to this problem is to find the constraints that these states satisfy, which set 
them apart from generic many-body states [3]. Such analysis can serve for the understanding of which type of 
entanglement that ground states can indeed harbour. To this aim, it is important to understand aspects of locality 
in these states. We ask: ‘to what extent can such states be described by a collection of local degrees of freedom, 
which are only loosely correlated with each other?’ 

Rigorous tools to tackle this question are scarce, even though various properties have been known in empirical 
ways (see [4]). An example is provided by the exponential decay of correlations, also known as exponential clustering: 
it has been proved that gapped ground states on a lattice have a finite correlation length, beyond which the 
correlations between spatially separated observables decay exponentially [5-7]. More recently, other quantitative 
tools have been devised, which characterize the ground state’s locality by looking at its entanglement structure [8, 9]. 
A notable example is area law of the entanglement entropy [8], which states that the entanglement entropy of a 
region with respect to the rest of the lattice should scale like the boundary area of the region rather than its volume. 
It is expected to hold for all gapped ground states on a lattice, but has only been rigorously proved in one spatial 
dimension (ID) by Hastings [10] (see Refs. [11-15] for further results). Hastings’ celebrated result yields a complete 
characterization of ID gapped ground states as matrix product states (MRS) [16], which, to a large extent, provides 
a full understanding of the ID case [17]. 

Unfortunately, in higher spatial dimensions our understanding of the problem is still very much limited. Not only 
that a proof for the area law is lacking, but it is also unclear how an area law would imply an efficient representation 
of the ground state [18[. Moreover, when the system has long-range interactions, or it is hosted in a lattice with 
a large dimensionality (like an expander graph [19]), locality properties of the ground state are even more illusive: 
exponential decay of correlations no longer holds (since all particles are essentially close to each other), and in 
general, area law become meaningless as surface areas become as large as volumes. For such systems, very well 
studied in the Hamiltonian complexity field, spatial distance might no longer a good figure of merit for identifying 
entanglement [20-22]. As we will shortly show, an alternative approach is to study entanglement and locality by 
analyzing the collective properties of a subsystem with respect to the number of local degrees of freedom it contains 
rather than the distance between them. 

In this paper, we introduce a new constraint on a many-body gapped ground states which complements some 
of the shortcomings of the existing approaches. We call it local reversibility. It is based on the intuition that 
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FIG. 1. Schematic picture of the local reversibility (LR). We disturb a quantum state |i/)) by an operator F^, which is 
supported in a subsystem L. We then try to recover the state Fn|'!/’) by the use of a g-local operator R. If the state |^/)) is a 
product state, we can recover the original state by an operator R with q = 0{y/\L\)\ then, ‘locally reversible state’ is defined 
as the class of states which have the same property as product state in terms of the non-locality of the reverse operator. 
The entanglement properties of LR states are expected to be highly restricted since entanglement cannot be recovered by 
local operations once it has been broken. 


macroscopic-scale entanglement cannot be recovered by any local operation once it has been broken. Therefore, 
states which allow this sort of local recovery, necessarily contain a ’small amount of macroscopic superposition’. 
Here, we observe that we use the term locality in a broader meaning than the usual spatial locality. 

We will show that such local reversibility holds for all unique gapped ground states of local Hamiltonians, includ¬ 
ing systems with long-range interactions or a diverging lattice dimensionality (for which the existing approaches to 
the locality properties, like the exponential decay of correlation, do not apply). We therefore believe that it exposes 
fundamental features of gapped ground states that cannot be captured by existing properties. To demonstrate 
its potential, we study specific problems in many-body physics. We work out rigorous bounds for the quantum 
fluctuations of locally reversible states. This, in turn, implies new constraints on the critical exponents and rigor¬ 
ous bounds on the quality of the mean-field ansatz, which is often used to treat complicated quantum many-body 
systems. An important outcome of our approach is an effective way to identify quantum macroscopic superposition. 


II. LOCAL REVERSIBILITY 


To motivate our approach, we begin with a heuristic discussion (Fig. I). Consider a state \ip) that is defined over 
N localized spins, each with a d-dimensional Hilbert space, and let F^ be an operator acting on a spin subset L; 
the total system is given by L U with the complement of L. Applying Yl to !■!/)), we can potentially disrupt 
the entanglement between L and L'’, even when Yl\iI}) has a constant overlap with |?/;). It is useful to think of \'tl)) 
as a superposition of several states !■!/)) = \'<pi) + \ 4 ’ 2 ) + ■ ■ ■ and of F^ as a projector that “kills” some (but not all) of 
these states. Intuitively, if l"^) contains some “global entanglement” on the scale of \L\ spins, we may only be able 
to reconstruct j-^) by acting on Y with (at least) an operator that acts non-trivially on the same portion L of 
the system (i.e., it would be an |L|-local operator). However, when |^) contains mostly short-range entanglement, 
we might be able to return to \'ip) by using an operator of a much smaller support. How much smaller should that 
support be for a slightly entangled state? Specifically, as we shall see shortly, the minimal size of support that is 
needed to reconstruct a product state is 0{^J\L\). This indicates that states that can be reversed by operators of 
0{\/\L\) support constitute a class of states with a small amount of entanglement. In the following, we refer to 
such a class as locally reversible states. 

We now put the discussion above on a formal ground. We first defines the notion of q-local operator, which may 
be often called a “few-body operator:” 

Definition II.1 (q-local) Given an integer q > 0, a q-local operator is an operator of the form O := X)|j 5 (:|<g 
where each ox is an operator supported on a finite subset of spins X = {ii, Z 2 ,..., i\x\\ of cardinality |A|. The ox 
operators are not necessarily sitting next to each other on the lattice. 

We formulate the reversibility property in terms of such operators ox- 

Definition II.2 (Local Reversibility) We say that a state \ip) is locally reversible (LR) if there exists a function 
f{x) that decays faster than any power law, such that for every subset of spins L and an operator Yl defined on it, 
and for every integer q > 0, there exists a q-local operator R such that 

WRTlW - |V’)II < |(J|r^|i/,)|/(^) ’ 


(1) 
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where || • ■ ■ || is the operator norm. 

Three remarks are in order, i) Both the shape and the size of L are left completely general. In particular, we 
can take L to be the entire system (|L| = N). ii) In some cases, it will make sense to only consider operators 
R that respect certain symmetries. We will later use this restricted definition of local reversibility for states with 
symmetry protected topological order (SPTO) [23]. The last remark is on the status of function f{x) in (1). 
Despite f{x) need to be a superpolynomially decaying function in II.2, the statement (1) itself can be proved for 
a fixed generic f{x). In this sense, the statement is non-asymptotic and valid for finite systems. In order (I) to 
be effective in putting bounds on the state in a meaningful way, however, f{x) need to be specific and non-trivial 
(see our main theorem III.l for an example of /(x)); such a feature will be thoroughly exploited in the rest of the 
paper. 

We claim that LR states show a specific degree of locality, while non-LR states correspond to states with non-local 
features due to global entanglement. This assertion can be explained by the following two lemmas characterizing 
the entanglement structure of LR states. 

The first lemma refers to the so-called macroscopicity of the states. Namely, we will demonstrate how non-LR 
states correspond to states with macroscopic superposition. 

Let’s consider states of the type |'0) = a\tjja) + /3|'*/’b) discuss the possibility that |'0a) and are macro- 
scopically distinct (meaning that a collection of local operators exists to IV’ft))- Then: 

Lemma II.3 Let |'0) be a state which satisfies (1) for a fixed function f{x). Then, for any decomposition \rp) = 
Prltp) + (I - Pl)\iP) ■= ctlfJa) + I3\ifb) with Pi = P^, we have 

OlV-a) = IV'b) + 1^), 

where O is a q-local operator and ||(5|p := (dj^) < f{q/^/\l7\). When \if) is a LR state, f{x) decays 

superpolynomially and only a difference of 0{ff'\L\) exists between the two states {tpa) o.nd {ifb}- 

The proof is provided in the Appendix A. 

By contraposition of the lemma, any quantum state such that we can find a bipartition P^ for which two states 
|'0 q) and li/jf,) are macroscopically distinct over a spatial scale, is non-LR; for example, the GHZ state over 
n particles, \tf) = ^(|0 • • • 0)„ -I- |1 • • • 1)„) is not LR since |0 • • • 0)„, |1 • • • 1)„ are clearly macroscopically distinct 
over the scale of n — 1, and we may write \if) = Pl\'4>) + (I — Pl)\'4’) with Pl := |0 • • • 0)(0 • • • 0|„. As we show later, 
this simple lemma also shows that degenerate topologically ordered states are not LR. 

The second lemma shows that fluctuations in an LR state are strongly suppressed. Indeed, consider an LR state 
\tf) together with a subset of spins L, and let be an additive operator of the form := Here, each Oi 

is an Hermitian operator with ||ai|| < 1, which acts only on the Ah spin. Since the Oi operators are commuting with 
each other, they can be viewed as classical random variables whose joint probability distribution is given by the 
underlying state The following lemma shows that their sum resembles a sum of independent random variables: 
its probability distribution is strongly concentrated around its mean with a width of 0{yJ\L\). 

Lemma II.4 Let and H^,,, be the projectors onto the eigenspaces of with eigenvalues < x and > x 

respectively, and let m be the median of A^ with respect to satisfying (1) in the sense that ('(/i|n<^|^) > 1/2 
and ('(/i|n>^|'0) > 1/2. Then, for any positive h the following inequality holds: 

(2) 

with a fixed function f{x). An equivalent statement is valid for ||. 

The proof is given by choosing Pl = in Lemma II.3. After a short algebra, we get ||n>^_|_^|^)|| < |/3| • 
\\^>m+hOAi.^\\+2f{q/y/\L\) with O g-local, where we use the facts |ap = > 1/2 and ||n^^_^^|V'h)|| = 

\\^>m+h\'^)\\l 1^- To finish the proof we will show that |ln>^_,_^On<^|| = 0 for q < h/2. This follows from the 
fact that A]^ is a sum of (commuting) 1-local operators of norm 1, and therefore every g-local operator can take 
an eigenvector |a) of A^ with eigenvalue a to a superposition of eigenvectors '^Cafa') with \a' — a| < 2q. Thus, 
choosing q = \h/2~\ — 1 proves the lemma. 

An immediate consequence of Lemma 11.4 with the assumption that fix) is a super polynomially decaying 
function, is that the fluctuations of every additive operator A^, which is defined on the entire systems (|L| = N) 
must satisfy 


((AA,.)^) := {^\Al\f:) - {fjlALW^ < OiN) . 


(3) 
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We point out that the well-known notion of macroscopicity measured by the Fisher information [24, 25] is implied by 
the Lemma II.3 and Lemma 11.4. This feature emerges clearly from the following reasoning. The Fisher information 
of a pure state |'0) with respect to an operator A is given by iFiip, A) = 4((AA)^) [25]. In Ref. [25] the authors suggest 
to define the ‘effective macroscopic size’ of a state as Wff(V’) •= max^ A)/(47V), where the maximization is 
over all extensive operators A := as in (3). States showing maximal quantum macroscopicity, such as the 

GHZ state, have Wff = 0{N), whereas states with no quantum macroscopicity have Nes = 0{1). Inequality (3) 
therefore implies that LR states have Wff = Equivalently, states with Wff = 0{NP) for p > 0 are necessarily 

non-LR. 

On the other hand, the converse is not true: there are states with Wff = £*(1) that are also non-LR. For instance, 
as we shall see, degenerate topologically ordered states turn out non-LR, but still satisfy the inequality (2), namely 
fVeff = 0{1). Thereby, LR provides us a more stringent characterization of the macroscopic superposition encoded 
in a many-body state. 


III. REVERSIBILITY OF GROUND STATES 


We now introduce our main tool for identifying LR states. The following theorem states that unique gapped ground 
states of local Hamiltonians are LR. It holds for a very wide class of quantum systems that are described by A:-local 
Hamiltonians of the form 

H = '^ hx with ^ \\hx\\<g Vi, (4) 

\X\<k X-.X^i 

where g is a constant of 0{1). Note that k is not necessarily equal to q from the definition of the operator R above. 
Also note that we implicitly assume that the spins sit on a lattice, but we make no direct use of the lattice structure 
or its dimensionality. Instead, we use the second condition in (4), meaning that the total strength of all interactions 
in which the Ah spin participates is bounded by a constant of 0(1). This definition of H captures a very wide 
class of quantum systems: with short-range interactions such as the the XY model, the Heisenberg model [26] and 
the AKLT model [27], as well as models with long-range interactions such as the Lipkin-Meshcov-Glick model [28] . 
Typically, we have k = 2 (i.e., two-body interaction), but several exceptions exist such as the ID cluster-Ising 
model [29] (fc = 3), the toric code model on a square lattice [30] (k = 4) and the string-net model on a honeycomb 
lattice [31] (k = 12). We denote the ground state of H by jH), and fix its energy to be Eq = 0. The rest of the 
energies are denoted hy 0 = Eq < Ei < E 2 < ■ ■ ■ ■ Finally, we let 6E := Ei — Eq be the spectral gap just above the 
ground state. With this notation at hand, our main theorem is given as follows. 


Theorem III.l With the above notations, for every spin subset L and every operator F ^ defined on it, and for 
any positive integer q, there exists a q-loeal operator R that satisfies 


where Uq := [q/k) and 

f ■= Ec = g\L\ + SgkuQ. 


(5) 

( 6 ) 


Inequality (5), together with the definitions of no and implies that IIRFi IH) — IH) II < and 

therefore |H) is LR when SE — 0(1). Hence the existence of a spectral gap places strong restrictions on structure 
of the ground states for very wide class of Hamiltonians. We note that the theorem requires no assumption on a 
spectral gap or the size of |L| and TV; hence, the theorem is not asymptotic and applicable for arbitrary ground 
states in finite systems. 

The full proof of Theorem III.l is given in Appendix B. Here we summarize its main ideas. Using recent results 
from Ref. [32], we conclude that after applying the operator F^ to the ground state IH), we get a state which 
consists mainly of excitations with energies of at most 0(\L\). Beyond that scale, the weight of the excitations 
decays exponentially. This is shown schematically by the blue curve in Fig. 2. Then following ideas from a recent 
new proof of the ID area law [12], we construct the operator R by approximating the ground-state projector using a 
polynomial of H. This polynomial is essentially a scaled version of the Ghebyshev polynomial (red curve in Fig. 2), 
chosen such that it approximately behaves as a boxcar function in the range [SE,2Ec + SE], thereby suppressing 
the majority of excitations in F^^IH). Grucially, even though it rapidly increases for x > 2Ec + SE, this blowup is 
cancelled by the exponential decay of the high-energy excitation. 
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FIG. 2. Schematic picture of the proof. After applying the operator to the ground state |f2), the energies at most of order 0{\L\) are 
excited (blue curve). We then filter out the excited states by an approximate boxcar function in the range [SE, 2Ec + SE] (red curve). 
Although the function rapidly increases for x > 2Ec + SE, this can be cancelled by the exponential decay of the energy excitation. 


Table 1: Locally Vs 

non-locally reversible states 

LR 

Non-LR 

Product state 

GHZ state 

Bounded-degree graph states 

States with large fluctuation 


Short-range entangled state Degenerate, topologically ordered ground states 

Degenerate, SPTO states (Symmetry-restricted non-LR) 


IV. EXAMPLES OF LOCALLY VS. NON-LOCALLY REVERSIBLE STATES 


Let us now apply Lemmas II.3, 11.4 and Theorem III.I to several exemplary states emerging in different contexts. 
The list of states is summarized in Table 1. In particular, we will demonstrate how local reversibility implies the 
absence of macroscopic superposition. We begin with LR states. 

1. Product states. A product state |'0) = IV'i) ® l'02) G) • • • G) IV'w) is LR because it is the unique ground state of 

the local Hamiltonian H = ~ \4’i){4’i \ ® Irest)- As H is made of commuting projectors, its spectral gap is 

necessarily 6E = 1. 

2. Graph states with bounded degree. These states are defined on a graph in which each node has at most 

0{\) neighboring nodes [33, 34] . The graph state is a non-degenerate gapped ground states of a Hamiltonian which 
is the summation of the following commuting stabilizers [35] {gi}^i: gt = erf (S> where [gi,gi'] = 0 

for {cr®, tJ^, cr^} are the Pauli matrices and {ji, j 2 , ■ ■ ■, jkt} are nodes which connect to the node i. By 

assumption, ki = 0(1), and hence the Hamiltonian is (!I(l)-local. By the commutativity of its terms, we conclude 
that it has a spectral gap SE = 0(1), and so by Theorem III.l such graph states are LR. 

3. Short-range entanglement (SRE) states. The third example are states that can be obtained by a constant- 
depth quantum circuit acting on a product state. In the literature they are often dubbed as “trivial states” [36, 37], 
or “short-range-entanglement (SRE) states” [2] . A constant-depth quantum circuit is a unitary operator that can 
be written as a product of fc = 0(1) unitary operators U = Ui ■ ■ ■ Uk where each unitary Ui is given as a product 
of unitary operators Ui = Ui^i ■ Ui ^2 • ■ • Ui^m with non-overlapping support of 0(1). To see why these are LR states, 
we write lip) = U\(j)), where U is the constant-depth circuit, and \(j)) = \(pi) ® [(^ 2 ) 0 • • • is a product state. Then 
it is easy to see that for any operator O with a support of 0(1), UOU~^ has also an 0(1) support, and therefore 
if iJ is a local Hamiltonian for which j^) is the unique ground state (see the first example), then El' = UE[U~^ 
is also a local Hamiltonian. Furthermore, El' has the same spectrum as H, and so it is gapped with the unique 
ground state, which is exactly \'ip). By Theorem III.l this state is LR. 

We note that not all LR states are also SRE states, or, equivalently, long-range entanglement (LRE) does not 
necessarily imply non-LR. For example, Kitaev’s toric code [30] on a sphere is a commuting local Hamiltonian and 
has a non-degenerate ground state with an 0(1) gap, and therefore by Theorem III.l it is LR. Nevertheless, it 
cannot be generated by a constant depth circuit working on a product state, and is therefore not an SRE state [38] . 
This point is also explained in Appendix C. 

We now turn to non-LR states. We will use Lemmas H.3 and H.4 to identify such states. 

5. “Schrodinger Cat” like states. States like the GHZ are not LR by Lemma II. 3. 

6. States with Fisher information of 0(Np) with p > 1. As we already mentioned, this result comes directly 
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from Lemma II. 4. Also here, a quintessential example of this class is the GHZ state [25], which has the scaling with 
p = 2. Moreover, the ground states at critical point are typically non-LR since they have p = 1 + {2 — rj — z)/D 
(see Appendix E), where 0 is the dynamical critical exponent, t] is the anomalous critical exponent, and D is the 
dimension of the system. For example, the critical point of the ID transverse Ising model has z = 1 and rj = 1/4, 
which yields p = 7/4. 

7. States with degenerate topological order. While the local fluctuations in Lemma II. 4 (as well as the Fisher 
information) cannot detect a locally hidden order such as the topological order, we can use Lemma 11.3 to see that 
states with a degenerate topological order are not LR. We demonstrate this point using Kitaev’s toric code model 
on a torus [30] with x ^/n sites. The idea is that by taking L to be a non-trivial loop in the torus of size y/n, 
there exists an operator Tl that takes one ground state |fli) to another ground state 1142), i-e., ID 2 ) = Tl\VI.i). The 
properties of the topological order guarantee that for any observable O that is supported on less than ^/n sites (the 
size of a Wilson loop), (r2i|0|r4i) = (r42|0|r42) and (r4i|0|r42) = 0. Therefore, we may invoke Lemma 11.3 with 
Pl := (I - Tl)/2, such that |Di) = Pl\^i) + (I - Pl)\^i) := a|r4+) + /3|D_), where |r4±) = ^(l^^i) ± ID 2 )). It is 

easy to verify from the above properties that 114+), |r4_) are macroscopically distinct over a scale of (they 

are in fact distinct over a scale of -^n, i.e., the size of L), and therefore by Lemma II.3, these degenerate ground 
states are not LR. 

We remark that non-degenerate topological order (e.g., in the toric code on the surface) results LR (from The¬ 
orem III.l). In this context, we observe that, despite the topological entropy is non vanishing for both degenerate 
and non degenerate topologically ordered ground states, the two cases are clearly distinct in terms of the irre¬ 
ducible multiparty correlation (the issue has been recently addressed in Refs [39-41]; see also Appendix C): Being 
our approach able to detect a ‘fine structure’ in the nature of the multipartite correlations, LR tells degenerate 
topological order apart from non-degenerate topological order. 

8. States with a degenerate symmetry protected topological order. The same arguments showing that 
degenerate topologically ordered states are not LR can be applied to the case of degenerate symmetry protected 
topological order. Such states show topological order only to a restricted set of operators defining a certain 
symmetry G [23]. They cannot be adiabatically connected to a product state using only operators from G, and 
in that restricted sense they are not SRE (see the following for the definition). An important example of states 
with SPTO can be obtained from graph state’s Hamiltonian on an open lattice, where one removes the boundary 
stabilizers. This removal introduces degeneracy to the groundspace. Much like the case of Kitaev’s toric code, we 
can also show here that the resulting ground states are non-LR as long as we restrict the operator R to satisfy 
the symmetry of the graph Hamiltonian without the boundary stabilizers. We refer to these states as symmetry 
restricted non-LR states. We present an example of such states for ID case [42] in the Appendix D. 


V. FLUCTUATIONS IN LOCALLY REVERSIBLE STATES 


Theorem III.l, together with Lemma II.4 provides a remarkable insight into the structure of unique ground 
statesL For any such ground state |f4), and for any additive operator A+ = defined on a spin subset 

I|n>m+;i|f4)|| < ^ with ci a constant of 0(1) (with m as defined in Lemma H.4). This implies that 

\{Al) — m\ = 0{\L\/6E), where (A+) = (f4|A+|f4) is the expectation of A+ in the ground state, and therefore 

l|n^(A,)+,|f4)|| (7) 

C 2 being a constant depending on the Hamiltonian’s parameters k and g. Taking A+ to be an order parameter (i.e., 
the magnetization in L), we arrive at the conclusion that the deviations of any order parameter from its expectation 
are exponentially suppressed in unique gapped ground states. It is interesting to contrast this inequality with the 
corresponding statistics of a product state. In such a case, A+ can be viewed as a sum of independent random 
variables, and by the Hoeffding’s inequality [43], lln>^^^^_|_^|'0)|| < /Rl). In this sense, unique gapped 

ground states enjoy a weaker, yet still non-trivial, notion of local independence. 

It is also worth noting that this independence cannot be (at least directly) deduced from the exponential decay 
of correlation of gapped ground states [6, 7], since it can be applied to sets of observables that may sit very 
close to each other on the lattice. Moreover, we can apply it to systems with long-range interactions, such as the 


^ We notice that because Theorem III.l is not asymptotic, the results in this section can be applied to arbitrary system size. 
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Lipkin-Meshcov-Glick model [28] and systems defined on the expander graphs [19], in which the maximal distance 
between any two spins is 0(1) and O(logA^), respectively. We remark that ineqnality (7) can be extended to 
generic few-body operators [44] : A = X)|x|<(} with q — 0(1); finally we can derive a similar bonnd for low-lying 
energy states, i.e., not necessarily the exact ground state (T.K., LA., L.A. and V. V., manuscript in preparation). 

A simple consequence of inequality (7) is a trade-off relationship between the spectral gap and the fluctuation 
AAl := ((f 2 |A||fl) - {niALinyy/"^ of a in the ground state: 

6E ■ {AAlY "A const • |L|,. ( 8 ) 


This has two interesting implications: 

1. Bounds on the critical exponents. As noted above. Theorem III.l does not assume the spectral gap of 0(1) 
and therefore can be applied to arbitrary ground states. Below, we apply it to quantum critical points to obtain a 
general inequality for critical exponents. 

Let us consider the critical regime, 6E —)■ 0. Define Aj^ = with L a total system and order 

parameters (e.g. magnetization). We then introduce the critical exponents z, rj, 7 and v as in Refs. [45[; z is the 
dynamical critical exponent, rj is the anomalous critical exponent, 7 is the susceptibility critical exponent and v is 
the correlation length exponent. By applying the finite-scaling ansatz [45] to ( 8 ), we can obtain 


. , 7 

Z > r — - = —, 
2 


(9) 


where the second equality comes from the Fisher equality 2 — rj = "f/v. We remark that (9) holds for very general 
settings both for homogeneous and disordered critical systems (see [46] for a non-trivial example where our inequality 
can be applied). Incidentally, we note that ( 8 ) gives non-trivial bounds for the critical Lipkin-Meshcov-Glick model, 
a system with long-range interactions [28, 47] . The details of this calculation are given in Appendix E 
2. Validity of mean-field approximations. Under the assumption of inequality ( 8 ) for ground states, we can 
estimate the validity of the mean-field approximation. Just as the first implication, the full details are given in 
Appendix F. The idea is that since the operators Al in ( 8 ) are arbitrary (as long as they are additive on L), we 
can use them to probe the two-spin reduced density matrix pij and its relation with its mean-field approximation 
Pi ® Pj- Specifically, it can be shown that for every spin subset L and an arbitrary spin i outside of it. 


\\pij - Pi® Pj II < const • ^J\L\/5E. (10) 

j&L 


This implies that on average, for each spin j € L, \\pij — pi ® Pj\\ < 0(1/^y\L\SE). If our system is defined by a 
nearest-neighbor two-body Hamiltonian on a regular grid with coordination number Z (the number of neighbors 
of each spin), then taking L to be the set of neighbors (|L| = Z), one immediately obtains a bound on the quality 
of the mean-field approximation for the energy density for VL 


1 

Z 


{hij)MF - 

<i,j> 


y~! {hij)e 


<i,3> 


< const • 


1 


where the sum is taken over the spins adjacent to i. We therefore obtain a quantitative bound on how the error 
of the mean-field approximation decreases as the lattice dimension (on which the coordination number depends) 
goes to infinity. This result is consistent with the folklore knowledge in condensed-matter physics that the mean- 
field becomes exact in infinite dimension. Recently, similar results have been obtained in different manners by 
Brandao et al. [48] and Osterloh et al. [49] In Ref. [48], the setup is more general (i.e., the system is not assumed 
to be gapped) but the error estimation is weaker than ours, scaling as 0{Z~^/^): In Ref. [49], the error estimation 
is as good as ours, 0 (Z“^/^), but under the additional assumptions of having a regular, isotropic, and bipartite 
lattice of |-spins. 


VI. SUMMARY AND OPEN QUESTIONS 

In this work, we introduced a new notion of locality in quantum states, the local reversibility, which is defined 
in terms of the type of local operations that are needed to reverse the action of perturbations to the state. 

We proved that all unique ground states of gapped local Hamiltonians are locally reversible (Theorem III.l), and, 
on the other hand, we showed how local reversibility implies a suppression of quantum fluctuations (Lemma 11.4). 
Together, these two results provide new insights into the structure of unique ground states of gapped local- 
Hamiltonians: i) a low Fisher information, which is an indication for the lack of quantum macroscopicity in these 
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states; ii) & novel inequality for the critical exponents in these systems; in) a quantitative analysis of the mean-field 
approximation; and finally, iv) since an adiabatic (local unitary) evolution of product states is locally reversible, 
our result clearly implies that all the gapped quantum phases of matter, disordered or with local order parameter 
(Landau symmetry breaking quantum phases), are reversible. In contrast, degenerate topological phases or the 
symmetry protected topological phases, are not reversible. We note that LR can detect the difference between 
degenerate and non degenerate topological order. Indeed, it was discovered that, although both with non vanishing 
topological entropy they have very different irreducible multipartite correlation (see paragraph 8 of SectIV and 
the Appendix C). In this context, we observe that LR can be further restricted (with a similar logic we pursued 
in this article to deal with symmetry protected topological phases) to improve and refine the characterisation of 
the ground state. Such a strategy might lead to catch properties of the state originating from the geometry of its 
ambient space. 

Our work provides an instrumental view for several research directions. 

Based on the bounds on the fluctuations we found, we might argue that, fluctuations in gapped ground state 
obey a Gaussian statistics (as they do in non interacting theories). A recent proof of the Berry-Esseen theorem for 
the quantum case by Brandao et al. [50] hints that this might be the case. A natural approach to this would be to 
tighten our main theorem, replacing the exponential decay in the RHS of inequality (5) by a Gaussian. 

Another intriguing direction to pursue is to incorporate LR, or one of its consequences, such as Lemma II.4 or 
inequality (8), - explicitly or implicitly - in the construction of tensor networks in higher dimension (e.g.. Projected 
entangled pair state, or PEPS [16]). By construction, these states satisfy the area-law, but we now know that they 
should also satisfy local reversibility. This will speed up the contraction of such tensor networks, which is the main 
bottleneck in the variational algorithms [51-54]. A goal of paramount importance in this context is to prove that 
PEPS are faithful representations of gapped ground states. A good place to start studying this question is in the 
ID world. We know that MPS can describe both LR and non-LR states (i.e., GHZ). The natural problem is then 
to pinpoint what is needed for an MPS to describe an LR state. 

Proving the area-law conjecture for gapped systems in 2D and beyond remains a challenge. It would be interesting 
to see if the additional structure of local reversibility of these states can assist in such proofs, or at least provide 
new insights regarding this important conjecture. As a specific route, we suggest to harness the LR in addition to 
the clustering, to improve the upper bound by Brandao and Horodecki [14] . 

Finally, it would be interesting to understand if local reversibility could somehow be used to characterize unique 
gapped ground states. In other words, is local reversibility also a sufficient condition for unique gapped ground 
states? Strictly speaking, this is incorrect, as there are LR states which are not gapped ground states. For example, 
the state jOOO • • • 0) -I- e(A^)|lll • • • 1) where e{N) decays faster than any polynomial is trivially LR, but can never 
be a unique gapped ground state of fc-local Hamiltonians as long as k < N/2 (see Ref. [55]). Nevertheless, we may 
still ask if, in some sense, every LR state can be approximated by a unique gapped ground state. If this is not the 
case, it would be interesting to understand which are these LR states that cannot be even approximated by gapped 
ground states. 

Generalising our approach to mixed states and devising experimental protocols to measure local reversibility are 
important future challenge. 


ACKNOWLEDGEMENT 

We are grateful to Naomichi Hatano, Tohru Koma, Hal Tasaki, Taku Matsui, Tomoyuki Morimae, Kohtaro Kato 
and Dorit Aharonov for helpful discussions and comments on related topics. We also thank Naomichi Hatano for 
valuable comments on the manuscript. This work was partially supported by the Program for Leading Graduate 
Schools (Frontiers of Mathematical Sciences and Physics, or FMSP) MEXT Japan, and World Premier International 
Research Genter Initiative (WPI) MEXT Japan. TK also acknowledges the support from JSPS grant no. 2611111. 
Research at the Gentre for Quantum Technologies is funded by the Singapore Ministry of Education and the 
National Research Foundation, also through the Tier 3 Grant random numbers from quantum processes. 


Appendix A: Proof of Lemma II.3 


Assume that fy) satisfies inequality (1). Then for every integer q > 0, there exists a g-local operator R such that 


RPlI'0 ) = W + where ||fy|p < 


IRi.ll 


^ Therefore, 


- V - lap 

/Jfyfc) = Ify) - RlW = (i? - I)PLfy) - l<5') 

= a{R-lMa)-\S'). 
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By denoting O = a{R — !)//? and |<5) = —j3 ^|(5'), we have IV'b) = 0\ipa) + |^), where ||5|p = ||^'||^/|/3| < 
completes the proof of Lemma II.3. 


Appendix B: Proof of Theorem III.l 
1. Outline 

The proof of Theorem III.l is rather technical, and therefore we first sketch it here, giving the full details in the 
following section. 

Multiplying inequality (5) by |(r2|rL|r2)|, and writing for brevity R := (0|rL|n)i?, we obtain 

||(i?-|0)(L!|).r,.|L!)|| <6||r,.||e-2"o/?. (Bi) 

So for the state to be LR, we need to find a R whose action on T/^jn) approximates the action of the ground state 
projector |n)(fl| on it. In addition, in order to satisfy the premise of the theorem, it has to be a q-local operator. 
To this aim, we look for a low-degree polynomial Ffi(x) and write R := Ffi{H). Specifically, choosing a polynomial 
of degree ng := \q/k\ guarantees that it will contain at most g-local terms, since, by definition, each term in F[ is 
fc-local. 

To understand the restrictions on Ffi{x) that inequality (Bl) poses, it is convenient to work in the energy 
basis expanding Tijn) = c(if)|if), we want i) Fji{0) = 1 (recall that have set Eq = 0), and ii) 

[J2e>se This is achieved using two ideas, which are demonstrated in Fig. 2. 

The first idea is that the expansion of Tijn) is dominated by energies of at most 0(|L|); beyond that scale, c{E) 
is exponentially decaying. This is a direct corollary of Theorem 2.1 in Ref. [32], which for our case implies: 

Corollary B.l (from Theorem 2.1 in Ref. [32]) Let he the projector into the eigenspace of F[ with ener¬ 
gies greater than or equal to E. Then 

Y, = ||nf^ri|fi)f < ||ri||2e-(®-25i^i)/4sfc. (B2) 

E'>E 

In Ref. [32] , this theorem was proved under the more restricted condition that every particle participates in at most 
g interactions of norm 1, but this can be easily relaxed to the current condition, given in definition (4). 

The bound in (B2) implies that our polynomial should mainly “kill” the energy excitations of Fijn) in the range 
\5E^O{\L\)\. Following Ref. [12], we let Eji{x) be the ngth order Chebyshev polynomial [56], scaled such that 
X : [—1,1] >—>■ [6E,2Ec -f 6E] and ^^■(0) = 1. As discussed in the following section, this polynomial fluctuates 
between in the range [SE,2Ec SE]^ and then diverges like 0{{2x/ Ec)’^°). It is our choice of Ec in 

Theorem III.l which guarantees that this divergence is cancelled by the exponential decay of Corollary B.l. After 
a rather straightforward calculation, one can show that total contributions of the energy segments \5E, 2Ec -f SE] 
and [2Ec + SE, oo) to || {R — |n) (f2|) • F || is exponentially small. 


2. Full proof 

Following the proof’s sketch in the previous section of the main text of the paper, we start from inequality (B2). 
Our goal is to find a polynomial Ffi{x) such that the action of the operator R Efi{H) on the state FijO) 
approximates the action of the ground state projector |0)(0| on it. As is a fc-local operator, choosing no := [g/fcj 
guarantees that ii is a g-local operator. 

Working in the eigenbasis of H, we expand Fl|0) = J^e '^iE)\E), and as Ffj(ii) is diagonal in this basis, 
[En{H)-\n){n\]TL\n) = {FR{o)-i)cio)\n)+ Y Fr{e)c{e)\e) . 

E>SE 

Therefore, for inequality (12) to hold, it is sufficient that 


Fr{0) = 1 

(B3) 



^ \c{E)Fr{E)\^\ <6||r,.||e-2"“/«. 

(B4) 

\e>5E / 
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As noted in the outline of the proof in the previous section, to prove these properties we use two ideas. The 
first is that the weight of the high energy excitations in rL|n) decays exponentially, as shown in Corollary B.l of 
Section B 1. The second is to take Ffi{x) to be a scaled version of the no’th order Chebyshev polynomial. Let us 
start from the second idea. The nth order Chebyshev polynomial [56] of the first kind is given by 


Tn{x) := 


(x + — 1)" + (x — \Jx^ — 1)" 

2 


(B5) 


Equivalently, for x G [—1,1] it is given by T„(a;) = cos(narccos(a;)), and for |a;| > 1 by Tn{x) = cosh(narccosh(x)). 
What makes the Chebyshev polynomial so useful to our purpose are the properties that are summarized in the 
following lemma, whose proof is given in Sec. B 2 a; 

Lemma B.2 


\Tn{x)\ < 1 , 

|r„w| < i(2|i|)” , 

\TM\ > 

Setting 


for |a;| < 1 
for |a;| > 1 

for |a;| > 1 





and Ec '■= g\L\ + SgkriQ , 


we define Fji{x) to be the polynomial 


Fr{x) 


rp ( X — SE 1 

Eg 


(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 


In other words, we defined it to be the noth order Chebyshev polynomial, scaled such that a: : [—1,1] i—>■ [5E, 2Ec + 
SE] and Fr{0) = 1. Clearly, this definition satisfies Eq. (B3). Let us see why it also satisfies inequality (B4). 

We begin by applying Lemma B.2 to the definition of Fr{x), which implies that for SE < x < 2Ec + SE, 


|Ffl(x)| < 2 e- 2 "«/«, 


(Bll) 


and for x > 2Ec + SE, 


\Fr{x)\ < 


2x - 2SE 


Er. 


- 2 


„- 2 no/C 


(B12) 


For brevity, we define the low and high energy ranges /low := \SE,2Ec + SE) and /hi := [2ifc + SE,oo). Then 
using the triangle inequality, we split the sum in the LHS of (B4) 

\ 1/2 / \ 1/2 / s 1/2 

^ |c(£;)Fi^(£;)n < E HE)FRiE)\A + E ’ 

^>6E / Vs'G/low / \E^Ihi / 


\E>SE 

and bound each term separately. The low-energy term is bounded by 

1/2 

2 e-" " ^ 


„-2«o/E |c(/;)| 2^ < 2e"^”“/«||rL|f2)|| < 2||rL||e"2”“/« , 

\£Je/LOW / 


(B13) 


which follows from Inequality (Bll) and the fact that Ee^Ii^ow < Ee = l|rL|fI)|P- 

To finish the proof, we will show that the high energies term is upper bounded by 4||rL||e“^"‘’/^. To this aim, 
we write /hi = /i U /2 U /a U ..., where Ij := [2Ec + SE + {j — \)g, 2Ec + SE -|- jg) and ry is a positive constant 
which will be set afterward. Using the triangle inequality once more, we get 

E |c(/;)Ffl(/;)p'j < E I E \ciE)FR{E)\^ 

EgIhi / j=i \EeIj 
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Clearly, for each Ij segment 


\ 1/2 

^ <max\FH{x)\ 

E^I, ) 


E I 


1/2 


As \Fji{x)\ monotonically increases for x > 2Ec + SE (which follows from the fact that the Chebyshev polynomial 
is monotonic for a: > 1 ), it follows that 


max|Fij(x)| < \Fr{2Ec + SE + jr]) \. 

xeij 


To bound the other term, we use Corollary B.l, which gives us 

1/2 / ^ 1/2 


,Eeij 


< 


iE>2Ec-\-SE-\-{j—l)r} 


where we have defined 


Together, this gives us 


A := 


Agk ■ 


1/2 


(B14) 


(B15) 


E HE)Fn{E)\^ < llrLlle^''- |Tj^(2A, + 5A + j77)|e-^(2E.+5£;+,r,-2g|L|)^ 

\E&i, J 

The final step is to show that for x > 2Ec + SE, 

|FR(a;)| • < g-2rao/? . ^-\ix-2g\L\)/2 

(see Subsection B 2b for a proof), which leads to 

\ 1/2 

E \c{E)Fii{E)\^ J < • e-^(2iJ,+5B+/r,-6ff|L|)/2 ^ 

Summing over all j > 1, then gives us 

/ \ 1/2 oo 

E HE)Fr{E)\'^ < ||rL||e"^"“/« • e-^(2i5,+5i5-2g|L|)/2 . ^Ar, ^ g-ar,A/2 ^ 

Vbg/hi / i=i 

Using the definition of Ec in Eq. (B9), we find that = e-i'(i6g*:"-o+'5B)/2 ^ calculating 

the geometrical sum we get E^i = e''‘’i/^/(l — which can be minimized to 4 by choosing r] 

such that e"’'''/^ = 2. All together, we therefore get 


^ \ 1/2 

E \c{E)Fr{E)\^\ <4||ri||e 

\Eeiiii / 


—2no/4 


(B16) 


which completes the proof. 


a. Proof of Lemma B.2 

Proof: Inequality (B 6 ) follows directly from the identity Tn{x) = cos(narccos(x)), which is valid for |x| < 1. For 
the other inequalities, first note that r„(—cc) = (—l)”T„(a;), which implies |T„(a;)| = |r„(|a:|)|, and so it is sufficient 
to prove inequalities (B7, B 8 ) for a; > 1. 
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To prove inequality (B7), consider the general inequality 

{2x - yr + y” < (2a;)" , (B17) 

which is valid for any a; > 1 and 0 < y < 1 (the inequality can be proved by differentiating (2a;)" — ((2a; — ?/)" + y") 
with respect to x, and noting for x > 1 and 0 < y < 1 it is a monotonically increasing function of x, and its 
minimum value 0, which is obtained for a; = 1 and y = 0). Choosing y = x — — 1, the LHS of inequality (B17) 

becomes 2T„(x), which proves (B7). 

For inequality (B8), we set t := arccosh(x), and then by the identity T„(x) = cosh(narccosh(x)), we conclude 
that for X > 1, 


r„(x) = cosh(nt) = i (e"* + e""*) > . 

To finish the proof, we need to show that for x > 1, t > 2w This follows from the fact that t/2 > tanh(t/2), 


and the trigonometric identity tanh(t/2) = ^ cosh(tj+ 7- 


b. Derivation of the inequality (B15) 


From inequality (B12), we have 


|^’i?(x)| < e-2"o/€ 


2x - 26E 


Er 


- 2 


for X > 2Ec + SE. To prove inequality (B15), we will show that [(2x — 26E)/Ec — 2]"°e 6g|i,|)/2 ^ ^ for 

X > 2Ec + 5E, or, equivalently, that its logarithm 


N A, „ ,,,, , /2x — 26E \ 

G{x) :=--(x- 65 |L|)+nolog(^- - - 2 j 


is negative. This follows from the facts that 


G[2E^ + 5E) = -2no 


\5E 


+ no log 2 < 0 , 


and for every x > 2Ec + 5E, 

dG{x) 

dx 


2 ^ x-Ec-5E - 2 ^ Ec 


A 

2 


A 

' no 


< 0 . 


Appendix C: Difference between degenerate and non-degenerate topological orders 

In the case of the toric code model, we find that the LR depends on the topology of the ambient manifold: LR 
holds on a sphere but is violated on non simply connected geometries (implying a non trivial ground-manifold). 
It is well-known, however, that the topological entanglement entropy is non-vanishing for toric code model ground 
states living in lattice with any topology [57, 58] . Indeed, the difference between the two kind of ground states can 
be resolved in terms of the irreducible multiparty correlation. 

The notion of irreducible multipartite correlation has been first introduced in Ref. [59] to characterize the mul¬ 
tipartite correlations in a quantum state. It was noted recently that such notion is equivalent to the topological 
entanglement entropy if the state has zero-correlation length [40] . As explained in Refs [39, 41], we have two 
kinds of multipartite correlation, which we refer to as ‘effective multiparty correlations’, distinct from ‘inherent 
multipartite correlations.’ The topological entanglement entropy cannot distinguish them. We have: 

i) The degenerate topological order, as that one of the toric code on a torus, has genuine multiparty correlation 
of the ‘inherent’ type involving 0{l) spins (1: the system length). 

ii) The non-degenerate topological order, as the toric code on the sphere, has low degree of inherent multiparty 
correlations involving 0{1) spins, but have the ‘effective’ type involving 0{l) spins. 
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^ _ System length : I ^ 



FIG. 3. Multipartite correlations in the surface code. In the ground states of the Kitaev model on sphere, it has no multi¬ 
party correlation (or contains only low-degree of correlations), but collective properties of the low-degree of correlations 
induce multi-party correlation when we look at reduced region of the system, say L'^. Indeed, if we split the region into 
A, B and C, we obtain non-trivial value of the topological entanglement entropy. 


In other words, a non-vanishing topological entanglement entropy in non-degenerate topological order arises just 
because of such multiparty low-degree correlations. There, we have no high-degree multiparty correlations if we 
look at the total system; in contrast, multiparty correlations of 0{l) can be effectively induced by tracing out some 
finite suregions (See Fig. 3) [41]. Such a conditional many-body correlations can appear in short-range entangled 
state [60, 61[ or even in classical models [62[. 

In this way, we can see qualitative difference between the degenerate and the non-degenerate topological orders 
in terms of the irreducible multiparty correlation, which results in LR of the surface code and non-LR of the toric 
code. Being our approach able to detect such ‘fine structure’ in the nature of the multipartite correlations, the LR 
tells degenerate topological order apart from non-degenerate topological order. 


Appendix D: Symmetry-restricted Local Reversibility 


Symmetry restricted LR states (SRL) can be introduced along very similar lines used in Sec. 11. Let’s consider a 
given Hamiltonian H enjoying a global symmetry G; let be the ground state of H. We say that the state ['(/') 
is SLR iff the property (1) holds with a g-local operator R enjoying the same symmetry group of the Hamiltonian: 
[R, G] = 0. 

Here we present an example of states which are not SLR. Cluster states provide an example of SPTO. The ID 
cluster states [63[ are the ground states of the Hamiltonian 




^i-1 


a, a 


X 

i+1 ’ 


(Dl) 


which enjoys a global symmetry Z 2 x Z 2 [42[. With the boundary conditions ctq = = I, the ground space 

of He is unique with a spectral gap. For ctq = cr£_|_i = 0, in contrast, the ground space is four-fold degenerate 
because the two stabilizers (out of L) crpcrjcr^ can be fixed at will [42[. Let {|n„),a = 0,1,2,3} 

be spanning the ground state manifold. Due to the symmetry-protected topological order of the system, it follows 
that the ground states |Dq) cannot be distinguished by any local operator ox in Z 2 x Z 2 : 

(flalojclDa) = (D^joxlfl/?), and (DqIoxIII/?) = 0. (D2) 

with |A| < cN (c = 0(1)). Using these conditions, the symmetry-restricted non-LR of jUo,) follows from the same 
arguments that were used in the proof of the non-LR of the toric code. 
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Appendix E: Critical exponents 

Here, we derive inequality (9) for the critical exponents z, 77 , 7 and v under the scaling ansatz (E2) [45, 64]. 
Recall that we are considering a local Hamiltonian system at T = 0 which is driven towards critically, and let 
A = ai, where ai are single particle operators that correspond to a local order parameter (e.g., spin localized 
at site i leading to the magnetization along a given axes). Our starting point is inequality ( 8 ), namely 

5E ■ (AAy < const • N. (El) 

We first define the variance (AHj)^ which depends on time as (AA^)^ := {{A{t) — (A)) ■ ((A — (A))), where 
A(t) = The variance (AA^)^ reduces to the summation of the correlation functions: 

JV N 

(AAt )2 = Mt)aj) - {ai){aj) := ^ Cij(t), 

where ai{t) := for i = l,2,...iV. Note that (AAt^o)^ is equal to (AA)^ = (A^) — (A)^. In the 

following, we denote Cij(t) = C{r,t) under the assumption of the translation symmetry. 

Now, we adopt the following scaling ansatz [45]: 

S{q,Uj;0 = e-^D{q^,ujn, (E2) 

where ^ is the correlation length and S{q,uj;^) is the spatial-temporal Fourier component of C{r,t), namely 

Siq,uj;0= j (E3) 

We also define S{q\^) as 

S{q-0 = ^J°^S{q,u;;0du^. (E4) 

We can see that the static fluctuation (AA(=o)^ is equal to NS{q = 0;^) by expanding S{q = 0;^). 

We then obtain the scaling of S{q = 0;^) oc i^y taking the scaling (E2) for S{q,uj-,^), and hence we 

have {AAi^qY/N oc We also have the scaling of the energy gap as 6Eq oc [45] by the use of the 

dynamical critical exponent z. At a critical point, where the correlation length is as large as the system length, the 
inequality (El) reduces to 

—z < —(2 — q — z) (E5) 

in the infinite volume limit (N —)■ 00 ). This reduces to the inequality (9) in the main manuscript. 

We close the section applying inequality ( 8 ) to a system with long-range interactions: the Lipkin-Meshcov-Glick 
model i?LMG = -^ with I 7 I < 1 . At the critical point A = \h\, we have the 

scaling [47] of SE oc and (AM^,)^ oc where is the magnetization in the x direction, = Y^=i ■ 

Thus, the spectral gap and the fluctuation can give the non-trivial sharp upper bounds to each other. 


Appendix F: The quality of the mean-held approximation 

Let jn) be the unique ground state of a gapped local Hamiltonian, and let pij,pi,pj be its two-particles and 
one-particles reduced density matrices. We want to estimate the error of the mean-field approximation pij —>■ pi®pj 
by proving inequality (10) in the main text. For simplicity, we set i = 1 and show that 

\\Pip - Pi® PjW < const ■ yJ\mdE . (El) 

First, note that we can always find a set of dP projectors onto the spin i = 1 that satisfy 

\\Pip - Pi® PjW < X] \\Pi^\pij - Pi® P])P[^'’\\ , 

m—1 


(F2) 
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where d is the local spin dimension. For example, in the case of spin-1/2 systems (d = 2), we can take = 
|0i)(0i|, = |li)(li|, Pf^ = |+i)(+i|, = |-i)(-i|, with |±i) := (|0i) ± \li))/^/2. Indeed, defining 

(5pij := Pi j - Pi ® pj, we get 

II^PijII < ll(0i|(5pip|0i)|| + ||(li|5pij|li)|| + ||(0i|(5pij|li) -I- (li|i5pij|0i)|| 

= ll(0i|(5pij|0i)|| + ||(li|^pij|li)|| + ||(+i|<^Pi,j|+i) - (-i|<5pij|-i)|| 

< ll(0i|(5pij|0i)|| + ||(li|^pij|li)|| + ||(+i|<^Pij|+i)|| + ll(-i|^Pi.i|-i)ll • 

The proof for higher d follows the same lines. 

Summing inequality (F2) over all j € L gives 

\\Pij - Pi<» PjW <YY - Pi Pj)Pi'^'’\\ ■ 

jGL m—lj£L 


To prove inequality (FI), we will show an upper bound of ~ Pi ® Pi)^i™^ll for arbitrary m. 

Defining p^™^ := Tri(P;^^™^pijP;^*'™^), where Tri(- • •) is the partial trace over the ith spin, we get 

- PI ® = pY^ ® (p^Y - ■ Pd) ■ 

Clearly, ||P;["*^ 0 (pj"*^ — • pj)\\ = \\p^Y ~ (f^lPi™^|f^) • Pj\\- Moreover, there always exists a rank-1 

projector such that 



. p,\\ = . Tr 


- > 




(D|p('")|fl).p,) 


where sign|Tr[pj"‘^(p^™^ — • pj)] |. Therefore, 


ll^i'"^(Pij - Pi ® Pj)^i'"^ll 
We now define the additive operator 


(m) rr 

= s; PTr 


(m) 
= ^1 


P^W(pW_(^|p(-)|^).p^.) 


(m) 




g(™) , pC 

j^L 


Then from the above calculation, 

Y W^Yipij - Pi ® Pj)P)'^^\\ = - (fl|Pi^'"^|D) • . 


But 


< ||P^”'^|fl)|l • - (D|y4(™)|D)|fl)|| 


and as |fi) — (D|y4(™^|n)|n)|j = AA^'^\ we conclude that 

Y \\Pi"^Hpi,j - Pi C) Pi)Pi^'"^|| < < const • y^\L\/SE. 

J&L 


Here, the last inequality comes from the inequality (8) in the main text, which applies in this case since is an 
additive operator on L. Combining this with inequality (F2) completes the proof. 
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Optimality of the bound 

When 5E = 0(1), inequality (FI) reduces to 

llPi.i - Pi ® PjW < const ■ \/\L\. (F3) 

J&L 

We can ensure that this upper bound is qualitatively optimal by considering the state 

^|0i)|0203 • • - 0^) + ^|li)|W2,..„^), (F4) 

where |W 2 ,.,,,Ar) is the W state for the spins 2,3,... ,N. We note that this state satisfies inequality (AA)^ < 0{\L\) 
[24], which is equivalent to the inequality ( 8 ) in the case of 6E = 0(1). Interestingly, the state in (F4) also gives 
the upper limit of the monogamy inequality of the entanglement [65]. 
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